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ABSTRACT 

We used Spitzer /IRS spectroscopic data on 426 galaxies including quasars, Seyferts, LINER and 
H II galaxies to investigate the relationship among the mid-IR emission lines. There is a tight linear 
correlation between the [Ne V] 14.3 /urn and 24.3 //m (97.1 eV) and the [O lV]25.9/iin (54.9eV) high- 
ionization emission lines. The correlation also holds for these high-ionization emission lines and the 
[Ne III] 15.56 /im (41 eV) emission line, although only for active galaxies. We used these correlations 
to calculate the [Ne III] excess due to star formation in Seyfert galaxies. We also estimated the 
O IV] luminosity due to star formation in active galaxies and determined that it dominates the 
O IV] emission only if the contribution of the active nucleus to the total luminosity is below 5%. We 
find that the AGN dominates the [O IV] emission in most Seyfert galaxies, whereas star-formation 
adequately explains the observed [O IV] emission in optically classified Hn galaxies. Finally we 
computed photoionization models to determine the physical conditions of the narrow line region 
where these high-ionization lines originate. The estimated ionization parameter range is -2.8 < log U 
< -2.5 and the total hydrogen column density range is 20 < log nn (cm -2 ) < 21. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: starburst — infrared: galaxies 



1. INTRODUCTION 

Primarily through measurements in the optical, ac- 
tive galactic nuclei (AGNs) have been categorized as 
Type 1 or 2, respectively with and wit hout very broad 
emission lines. The un ification model (|Antonuccil 119931 : 
lUrrv fc Padovanilll995|) successfully explains this behav- 
ior in terms of a dusty circumnuclear torus that hides 
the broad-line region (BLR) from our line of sight for 
type 2 AGN. It also appears th at extinction in the host 
galaxy can hide the BLR (e.g., iMaiolino fc RleklU995t 
lAlonso-Herrero et~a l. 2003) in some active galaxies. Ex- 
tinction therefore affects the optical properties of AGN 
in fundamental ways. These effects can be minimized by 
studying these objects using infrared (IR) emission lines, 
both to test the predictions of the unification model and 
to characterize AGNs in a uniform way Toward this 
end, we have used mid-IR observations with the Infrared 
Spectrograph (IRS. lHouck et al.ll2004T ) on Spitzer to ex- 
amine the behavior of a large sample of AGN. 

Important high-ionization lines accessible in the mid- 
IR include [Ne V] (97.1 eV) at 14.32 and 24.32 /mi, 
[O IV] (54.9 eV) at 25.89 /xm, and [Ne III] (41 eV) at 
15.56 /im. Because of their very high ionization po- 
tential, the [Ne V] lines are considered to be reliable 
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signp osts for an AGN (|Genzel et al.l 119981 : lArmus et al.l 
2007). These lines have been used to estimate the ac- 
cretio n power in the local Universe ([Tommasi n et al.l 
12010]) and to identify low- l uminosity AGNs in local galax - 
ies ([Satvapal et al.l [20(51 iGoulding fc Alexand"en l2009( ) . 
Howev er, they are also p r oduced in supernov a remnants 
(SNR, IQliva et al.1 [19991: ISmith et all 120091) . planetary 



nebul ae (PN,| Pottasch et al 



(WR, iSchaerer fc Stasiriskal 



2009f) . and Wolf-Rayet stars 



1999). Although the [Ne V] 



luminosities of these objects are low, ~ 10 erg s 
(jSmith et al.l 120091: iPottasch et~aTI I2009D . several thou- 
sands of them might also produce detectable [Ne V] emis- 
sion. 

Despite the lower ionization potential of [O IV], the 
25.89 /im emission line appears to be an accurate in- 
dicator of AGN power, since it corre lates well with 
the hard (> 14keV ) X-ray luminosity dMelendez et all 
I2008ai iRigbv et al.l 120091: iDiamond-Stanic et al.l 12009ft 
and t he mid-IR [N e V] e mission lines (jDudik et al.l 
120091 : iWeaver et all 120101 ). However, this line 
also appears in the spectra of starburst galaxies 
with no other evidence for AGN (jLutz et alj Il998at 
iBernard-Salas et al.ll2009|) . where i t is attributed to Wolf- 
Raye t stars (jCrowther e t al. 1999; Schae rer fc Stasihska 
1 19991) and/or to shocks (jAUen et al.1 120081 : ILutz et al. 
Il998al) . 

Th e [Ne HI] 15.56 /im line is excited by young, massive 
stars dVerma et al.ll 2003: Br andl et al . 2006: Beira o et al 



2006, 20 081: iHo fc Ketol l2007t [B ernard- S ala s et al 



2009; Alonso-Herrero et alj 



2007; 

mm 
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120101 ). Nonetheless, it seems to be a reasonably good 
proxy for AGN luminosity, a t least for reaso nably 
high luminosities dGorjian et al.l[2007l iDudik et al.ll2009l : 
iMelendez et al.ll2008al : iTommasin et al.ll2008j ). Quantifv- 
ing when its excitation is dominated by an active nucleus 
would allow it to be used in concert with the other high- 
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ionization lines to probe conditions in the neighborhoods 
of AGNs. 

In this paper, we present a statistical study of the be- 
havior of these lines in a sample of 426 galaxies. The 
data are compiled from the literature on high spectral 
resolution (i? ~ 600) IRS measurements, as well as our 
own new measurements of low-luminosity Seyfert galax- 
ies and luminous infrared galaxies (LIRGs). We describe 
the sample and the data reduction respectively in Sec- 
tions 2 and 3. Section 4 explores the correlations among 
the high-ionization emission lines. In Section 5, we dis- 
cuss the star formation contribution to these lines. Sec- 
tion 6 compares the star formation in Seyfert 1 and 2 
galaxies. Finally, in Section 7 we use photoionization 
models to study the physical conditions in the narrow 
line regions (NLRs). 

Throughout this paper we assume a flat cosmology 
with H = 70 km s^Mpcr 1 , fi M = 0.3, and fl A = 0.7. ' 
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Fig. 1. — Distribution of the galaxy distances of the complete 
sample (black), QSO (gray), Seyfert 1 and 2 galaxies (red), LINERs 
(blue) and Hll galaxies (orange). 



TABLE 2 
The sample 



2. THE SAMPLE 

The sample contains 426 galaxies (Table [lj for which 
high spectral resolution (R ^600) Spitzer /I RS spectra 
were available, either i n the literature dWeedman et alJ 
20051: lOde et all 120061: iFarrah et all 120071 : iGoriian et al. 
2007 [ [Tommasin et all I2008L 120101 IVeilleux et al l 
200a [Dale et alJ 120091 iBernard-Salas et alJ 120091: 
Goulding k Alexander! (2001 Pereira-Santaella et al. 
2011, in prep.), in the Spitzer archive or observed 
through the programs 40936 and 50597 (PI: G. H. 
Rieke). AGN type classifications were taken from 
NASA Extragalactic Database (NED). We designated 
only pure cases as Type 2 and included within Type 
1 all intermediate cases (i.e.: 1.2, 1.5, 1.8, 1.9) based 
on observations and mod eling of the nuclear spec- 
tra energy distributions rtAlonso-Herrero et al.l 120031 : 
I Ramos Almeida et alJ I2009I ). We could not find the 
nuclear activity classification for 32 of the galaxies. 
The sample includes 28 QSOs, 76 Seyfert 1, 125 
Seyfert 2, 55 LINERs and 110 Hll or starburst galaxies 
(see Table [2]). The luminosity of the active galaxies 
range s from QSO (QUEST sample, IVeilleux et afl 
20091) to t ypica l Seyfert galaxie s (12 yin sample, 



Rush et aLI Il993t iTommasin et al.l 
low-l u minosity Seyfer t s (RSA sample, 



2010 ), relatively 



Maiol 



119951: IHo et all Il997t iDiamond-Stanic et al.l 120091 ) and 



no k Riekd 



LINERs (|Sturm et al.H2006f) . According to their infrared 
luminosities 71 sample members are classified as LIRGs 
(-Lib, = 10 11 — 10 12 L Q ) and 54 as ultraluminous infrared 
galaxies (ULIRGs, L m > 1O 12 L ). We found in the 
Revised Bright Galaxy Sample catalog (| Sanders et alJ 
12001 : ISurace et al.ll20f)l the IRAS fluxes for 196 of our 
ga laxies which we used to ca lculate their Lir as defined 
in lSanders k Mirabell (|1996t ). 

Note that the LINER group includes very different 
galaxies. IR-bright LINERs (Lir/Lb ^ 1) have infrared 
SEDs similar to starbursts, although high-ionization lines 
and hard X-ray cores are detected in some, suggesting 
the presence of an AGN. In comparison, IR-faint LIN- 
ERs (Ltr/Lr < 1) seem to be power ed by an AGN 
(jSturm et alJl2006t ISatvapal et al.ll2004D . 

Figure [1] shows the distance distribution of the sample. 
The median distance is 60 Mpc and most (68%) of the 
galaxies are between 15 and 400 Mpc. 



Type 


N 


[Ne V]14.3^m 


[Ne V]24.3^m 


[O IV] 


QSO 


28 


22 


16 


25 


Seyfert 1 


76 


59 


42 


63 


Seyfert 2 


125 


89 


69 


101 


LINER 


55 


10 


6 


27 


Hn/Starburst 


110 


4 


2 


59 


Unknown 


32 


2 


1 


7 


Total 


426 


186 


136 


282 



Note. — N is the number of galaxies of each type. For each 
type we give the number of detections of the [Ne V] 14.32 fim, 
[Ne V]24.32 /un and [O IV]25.89 /an emission lines. 

3. DATA REDUCTION 

In addition to results from the literature, we utilized 
the measurements for the 88 galaxies of the RSA sam- 
ple observed with the high spectr al resolution modules 
as wel l as for the 34 LIRGs of the lAlonso-Herrero et alJ 
(2006) sample (3 of these LIRGs are also members of 
the RSA sample) . We retrieved the basic calibrated data 
(BCD) from the Spitzer archive processed by the pipeline 
version S18.7. We subtracted the background contri- 
bution when a dedicated sky observation was available. 
Note however that the sky subtraction is not important 
to measure these fine structure emission lines and it just 
improves the final quality of the spectra by removing bad 
pixels. Then we extracted the spectra using the standard 
programs included in the Spitzer IRS Custom Extraction 
(SPICE) package provided by the Spitzer Science Center 
(SSC). We assumed the point source calibration for all 
the galaxies, which is a good approximation for most of 
the galaxies at least for the high ionization lines. For 
the galaxies observed in the mapping mode we extracted 
the nuclear spectra from the data cubes using a square 
aperture of 13'.'4x 13'.' 4 and then we applied an aperture 
correction (Pereira-Santaella et al. 2011, in prep). We 
used a Gaussian profile to fit the emission lines. For 9 
galaxies (NGC 777, NGC 3254, NGC 3486, NGC 3941, 
NGC 4138, NGC 4378, NGC 4472, NGC 4698, 
NGC 5631) we were not able to measure any spectral 
feature due to the low signal-to-noise ratio of the spec- 
tra. They are not listed in Table [1] nor included in the 
analyzed sample. 

Table [2] includes the number of detections of the 
[NeV] 14.32 /jm, [Ne VJ24.32 fim and [O IV]25.89/mi 
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TABLE 3 
Median line ratios 







Median ratio 




Ratio 


QSO 


Syl Sy2 LINER 


Hii 



O IV] 25.89 /im/[Ne V] 24.32 fim 3.3±1.2* 3.6±0.6 3.5±0.6 3.7±1.0* 4.1±1.3* 

O IV]25. 89 /ttm/fNeV] 14.32 fim 3.6±0.9 3.5±1.0 3.5±1.3 4.8±2.0* 4.4±2.0* 

Ne V]24.32^m/[Ne V] 14.32 /an 1.0±0.3* 1.1±0.3 1.0±0.3 1.3±0.3* 4.7±4.1* 

[Ne III] 15.56 Aim/ [Ne V]14.32^m 1.7±0.5 1.9±0.6 2.0±0.7 5.3±3.7* 6.8±1.7* 

[NeII]12.81/xm/[Ne V]14.32Aim 0.7±0.4 1.6±0.7 2.0±1.1 28±14* 57±42* 

[NeIII]15.56Aim/[OlV]25.89Atm 0.5±0.2 0.6±0.2 0.9±0.4 2.3±1.1 4.2±2.4 

[NeII]12.81/im/[OlV]25.89Atm 0.2±0.1 0.5±0.3 1.1±0.9 6±2 22±16 

[NeIII]15.56Atm/[Ne II] 12.81 Aim 2.0±1.0 1.2±0.6 0.8±0.4 0.3±0.2 0.17±0.07 



Note. — Median ratios and uncertainties for each type of galaxies. The uncertainty 
is calculated as the median absolute deviation (should be multiplied by 1.48 to obtain 
the standard deviation). 

* These values are calculated with less than 20 galaxies. 



emission lines and in Table |3] we show the observed me- 
dian line ratios for each galaxy type both for galaxies 
from literature and those analyzed by us. Note that 
there are 4 detections of the [Ne V] 14.32 Aim line and 2 
of the [Ne V]24.32/im line in galaxies classified as Hn. 
For these galaxies (NGC 613, NGC 1792, NGC 3621 and 
NGC 5734) the detection of the [Ne V] lines is the only 
evidence of AGN activity. We decided to keep their H n 
classification because: (1) the AGN may be extremely 
obscured or very low luminosity and, thus, the nuclear 
spectra might be dominated by star formation features; 
and (2) these lines are also detected in SNR, PN and WR 
stars which could be the origin of the [Ne V] emission in 
these galaxies. 

Table |3] gives the median and deviation of all the line 
ratios we discuss in this paper for each type of galaxy. 
The line fluxes for all the galaxies are listed in Table [TJ 

4. THE MID-IR HIGH-IONIZATION EMISSION LINES 
4.1. The [0 IV] versus [Ne V] correlations 

The [Ne V]24.32Aim and the [O IV]25.89 Aim emis- 
sion lines are commonly detected in active galaxies 
(iLutz et al.lll998ri iGenzel et all 119981: ITommasin et al.l 
|2010() . The detection of the former is considered a good 
indicator of AGN activity because of its high ioniza- 
tion potential. However using the [O IV]25.89 /an line 
as an AGN tracer is not as straightforward. For opti- 
cally classified Seyfert galaxies the [O IV]25.89/im lumi- 
nosity is a good proxy for the AGN intrinsic lumin osity 
(|Diamond-Stanic et al.|[2009t IMelendez et~aT1l2008aD and 
there is a good correlation between the [N e V]24.32 Aim 
and t he [O iV]25.89Aim e mission in AGNs (jDudik et al.l 
[2(5ol iWeaver et al.l[20To[ ). Howev er, the [O IV] 25.89 Aim 
line i s also produced by WR stars (jSchaerer fc Stasinskal 
119991) and is observed in starburst gal axies with no other 
evidence for the presence of an AGN (ILutz et al.lll998al 
iBernard-Salas et al.|[2009l ). 

The left panel of Figure [2] shows the tight correlation 
between the [Ne V]24.32 Aim and [O IV]25.89Aim lines 
spanning at least 5 orders of magnitude in luminosity 
(38 < log i[o iv] ( er g s_1 ) < 43). We do not find any de- 
pendence between the optical classification of the nuclear 
activity and this correlation. This implies that the ion- 
ization parameter is similar in type 1 and type 2 Seyfert 



galaxies (see Section [7]). The slope of the best fit to 
the luminosity data is 0.96 ± 0.02 with a 0.14 dex dis- 
persion. It is reasonable to assume a linear correlation 
between these two emission lines, and this yields a ratio 
[O IV]25.89 Aim/[Ne V]24.32 Aim =3.5, with rms scatter of 
0.8. 

The [Ne V]24.32Aim line is detected in -50% of the 
optically classified Seyfert galaxies (see Table [2]). How- 
ever, for a considerable number of galaxies we have the 
[O IV]25.89Aim measurement and an upper limit for the 
[Ne V] 24.32 Aim flux. We plot these fluxes and upper lim- 
its (middle panel of Figure [2]) for all the galaxies. The up- 
per limits to the [Ne V] 24.32 Aim flux are compatible with 
the correlation for most of the Seyfert galaxies. However, 
the [O IV]25.89Aim line is detected in more than 50% of 
the H II galaxies in our sample. In 90% of these galaxies, 
the sensitivity of the [Ne V]24.32 /im measurement (or 
upper limits) is inadequate to probe whether the [O IV] 
emission is associated with a hidden AGN. 

In short, we find that the [O IV]25.89/im luminosity is 
well correlated with the intrinsic AGN luminosity (i.e., 
[Ne V] luminosity) for Seyfert galaxies and quasars with 
L[o iv] > 10 39 erg s _1 to, at least, L[ iv] ~ 10 43 erg s" 1 . 
For lower luminosities the fraction of the [O IV]25.89 /im 
emission produced by star formation may be consid- 
erable. Thus the [O IV]25.89 /im may not be an ac- 
curate tracer of the AGN luminosity for galaxies with 
L [Q Iv] < 10 39 erg s- ^Goulding fc Alexander] l2l)09l ). 

The linear correlation also holds for the [Ne V] 14.32 Aim 
and the [O IV]25.89/im luminosities (Figure[3]), although 
the dispersion is larger ([O IV]25.89 Aim/[Ne V] 14.32 Aim 
=3.4, , with rms scatter of 1.4 assuming a linear fit). 

4.2. The [Ne III] and [Ne II] emission lines 

Due to the intermediate ionization potential of 
Ne+ 2 (41.0eV), the [Ne III] 15.56 Aim emission line 
may be produced by young stars or by AGNs. 
This line i s obs e rved in sta r form i ng galaxies 
dVerma et all 120031: iBrandl et al.l 120061: iBeirao et all 
20061 120081 : IHo fc Ketdl2007t IBernard-Salas et all 120091: 
Dale et all l2009t iPereira-Santaella et all I2010D , and it 
is also correlated w i th the AGN in t rinsic luminosity 
dGoriian et all 120071: iDeo et al.l 120071: ITommasin et all 
20081: IMelendez et alll2008bD . 




Fig. 2. — Relationship between the [O IV]25.89^m and [Ne V]24.32 /im luminosities (left) and fluxes (middle). The red line is the best fit and the black line the best linear fit. The 
dashed lines mark the lcr deviation. QSO are plotted as gray circles, Syl as red circles, Sy2 as green circles, LINERs as blue circles, Ho galaxies as orange circles and thos e galaxies 
without a classification as purple circles. The open symbols mark the galaxies with star formation which may contribute to the [O IV]25.89/xm emission (see Section 14.21 and right 
panel of Figure|4]|. The right panel shows the [O IV]25.89 /Ltm/[Ne V]24.32^tm ratio vs. the [O IV]25.89/^m luminosity. The solid red line is the ratio obtained from the linear fit and 
the dashed lines mark the lcr deviation. 




Fig. 3. — Relationship between the [O IV] 25.89 fim and [Ne V] 14.32 fim luminosities (left) and fluxes (middle). The right panel shows the [O IV] 25. 89 /im/[Ne V] 14.32 (im ratio vs. 
the [O IV]25.89/xm luminosity. Symbols are as in Figure[2] 




Fig. 4. — Relation between the [Ne III] 15.56 fim and [Ne V] 14.32 /an emission. Galaxy symbols are as in Figure[2] The open symbols mark those galaxies above 3<r the [Ne V] 14.32 /an 
vs. [Ne III] 15.56 fim flux correlation. The black line in the left panel is the ratio obtained from the flux correlation. In the left panel only are included the [Ne V]24.32 fim upper limits 
of the QSO and Seyfert galaxies. 
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Fig. 5.— [Ne Il]12.81/im/[Ne V]14.32/im ratio 
[Ne V] 14.32 /jm luminosity. Symbols are as in Figure[4] 



The middle panel of Figure 0] shows the relation be- 
tween the observed flux of the [Ne III] 15.56 /im and 
[Ne V] 14.32 /im emission lines for our sample of AGNs 
and for the 4 Hn galaxies with [Ne V] detections. We 
compare these lines because they are close in wave- 
length and the effect of differential extinction is min- 
imized. It is clear that some galaxies have an ex- 
cess of [Ne III] 15.56 /im emission relative to that of the 
[Ne V] 14.32 //m, which we attribute to star formation (see 
Section [5]) . Therefore to obtain the fit to the data we 
use the outliers resistant linear fit algorithm provided 
by the IDL function ROBUST.LINEFIT. The slope of 
this fit is 1.02 ±0.02 with a scatter of about 0.2 dex 
around the fit. The result is consistent with a simple 
correlation; by fixing the slope to unity, we obtain the 
ratio [NeIII]15.56Atm/[Ne V]14.32 / um =1.9, , with rms 
sca tter of 0.8. This corr elation was previously reported 
by IGoriian et al.l (|2007t ) , although they found a slight 
dependence of the ratio with the [Ne V] 14.32 /im lumi- 
nosity (slope 0.89). Their sample included 53 X-ray se- 
lected AGNs with [Ne V] 14.32 /im luminosities between 
10 39 and 10 43 erg s _1 . The larger number of galaxies 
in our sample allows us to minimize the contribution of 
the galaxies contaminated by star formation to the fit 
and the luminosity dependence of the correlation is sub- 
stantially reduced. Where the AGN clearly dominates, 
the points fall around the line determined for the whole 
sample by ROBUST_LINEFIT, thus supporting our in- 
terpretation that this fit determines the relationship for 
AGNs in general. 

The [Ne II] 12. 81 um emission traces young (<10 Myr) 
star formatio n dRoche et al.1 119911: iThornlev et al.l l2000t 
Verma et all 120031: iRigby fc Riekei l2004t iSniiders et al.1 
20071 IHo fc Ketol l2007t iDiaz-Santos et al.1 12010ft . The 



ionization potential of this emission line is 21 eV and thus 
it is mainly produced in H n regions. Figure [5] shows 
that the scatter in the [Ne II]12.81 //m/[Ne V]14.32/im 
ratio (0.7 dex) is larger than that found in the other 
line ratios, for the lower AGN luminosities there is a 
large scatter. There is a weak trend of [Ne II] 12.81 /im/ 
[Ne V] 14.32 /im ratios decreasing with increasing AGN 
luminosity. Assuming that there is an intrinsic ratio 
for AGN, this behavior implies that the AGN does not 
dominate the [Ne II] 12.81 fim output and that there is 
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Fig. 6.— [O IV]25.89A*m/[Ne Il]12.8lAim ratio vs. 
JO IV]25.89 fim luminosity. Galaxy symbols are as in Figure 
[2] The dashed orange ([O IV] /[Ne II] = 0.05) and green ([O IV]/ 

[Ne II] = 0.35) lines are the boundaries of the Hn and Seyfert 

galaxies respectively. 

a wide range of star formation rates contributing to 
its luminosity. However, for the most luminous objects 
(logL[ Ne vji^m (erg s _1 ) > 41.5), the smaller scatter in- 
dicates that the relative contribution of the AGN to the 
t otal [Ne II] 12. 8 1 ixm excitation is higher. 

iSturm et all (|2002f) determined that the pure AGN 
[Ne V] 14.32 /xm/ [Nell] 12.81 /Ltm ratio is 1.1, but we find 
many Seyfert galaxies with higher ratios. The largest 
[Ne V]14.32/im/[Ne II] 12.81 /im ratio in our sample is 
~5.5, which is compatible with the larg est ratio pre- 
dicte d by NLR photoionization models (jGroves et al.l 
120061 their Figure 11). With our data (Figure [5]) we 
are not able to determine the pure AGN ratio, which 
also d epends on the ionization parameter (jGroves et al.l 
I 2006D. but it is lik ely that it is between that found by 
ISturm et al.l (|2002l) and the largest ratio in our sample. 

The [O IV]25.89/xm/[Ne II] 12.81 ^m ratio can be 
used to separate A GNs from star-forming galaxies 
dGenzel et alj 119981: ISturm et al.l l2002t iPeeters et all 
120041: iDale et al.l I2009Q . We plot this ratio versus the 
[O IV]25.89 /im luminosity in Figure El Most of the 
active galaxies have ratios larger than 0.35 whereas 
the ratios for Hu galaxies are lower than 0.05 (Fig- 
ure 5 of IDale et al.l 120091) . LINERs appear in the 
region between these ratios (Figure [6|. A consider- 
able number of Hn galaxies also have [O lV]25.89/im/ 
[Ne II] 12.81 /xm ratios above 0.05. The [Ne V] lines are 
not detected for them although their upper limits are 
compatible with the [O lV]25.89/mi versus [Ne V] cor- 
relations and thus a low-luminosity AGN could be the 
origin of the [O IV]25.89/im emission. The large scat- 
ter in the [O IV]25.89/im/[Ne II]12.81yiim ratio, due to 
the different star formation rates in the AGNs, does not 
allow us to determine the value of the pure AGN ra- 
tio. The largest [O IV]25.89/im/[Ne II]12.81 /im ratios 
that we find in AGNs are ^10. 

5. STAR FORMATION CONTRIBUTIONS TO THE MID-IR 
HIGH-IONIZATION EMISSION LINES 

5.1. [O IV] Contamination by Star Formation 

In Figures O and [3] there are few outliers above the 
[O lV]25.89/im versus [Ne V] correlations. Most of them 
are represented by open symbols, which indicate that 
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their [Ne III] 15.56 /im/[Ne V] 14.32 fim ratios (see Sec- 
tion 15. 2p are larger than those found in AGN domi- 
nated galaxies (QSO, Syl and Sy2). It is likely that 
the [O IV] 25.89 /im emission of these galaxies is con- 
taminated by star formation. We now compare the 
[O IV] 25.89 /im emission of star forming regions and 
AGN. The [O IV]25.89 /xm/[Ne II] 12 81 fim ratio in star 
forming regions is < 0.05 (Section l4~2l . Taking this upper 
limit and using the relation between the [Ne II] 12.81 fim 
luminosity and the tota l infrar ed (8-1000 /im) luminosity 
(L m ) from IHo fc Ketol (pOOl we obtain the following 
relation for star-forming galaxies. 

logL [0 IV] (erg s- 1 ) < logi m (erg s^ 1 ) -4.7 ±0.6 (1) 

Rigb v et"al~l ()2009P ) calculated the ratio between the 
total AGN luminosity and the [O IV] 25. 89 /im emission. 
We use the ratio for type 1 AGNs (~2500) because the 
hard X-ray (E > 10 keV) emission used for the bolometric 
corrections may b e affected by extinction in Seyfert 2s 
(|Rigbv et al.lf2009h . 

\ogL lo IV] (erg s" 1 ) = \ogL AGN (erg S - 1 )~3A±0A (2) 

Then, combining both relations we find that the lumi- 
nosity due to star formation, Lir, has to be at least 20 
times brighter than the AGN luminosity (that is, AGN 
contribution to the total luminosity below 5%) for the 
star formation to dominate the [O IV]25.89/im emission. 

We note that the number of H II galaxies falls rapidly 
above log L[o i V ] ( er g s_1 ) =40.2 (see Figure |6|). Assum- 
ing that star formation dominates the [O IV]25.89/im 
emission at this luminosity in these galaxies, the equiva- 
lent infrared luminosity (Lir) would be 2xlO n L©. This 
value of the Lri seems reasonable since most of the H n 
galaxies with log L[ iv] ( er g s_1 ) between 39.6 and 41.0 
are classified as LIRGs. The Lir limit c alculated above 
is one order of magnitude larger than L± (jTakeuchi et al.1 
2003) so the reason for the lower number of H II galax- 
ies above this [O IV]25.89 /im luminosity can be that the 
space density of luminous star forming galaxies decreases 
rapidly with increasing luminosity. 

Out of 196 galaxies with estimates for Lir, 
[O IV]25.89/mi is detected in 122. Only 2 ±2% of the 
Hll galaxies have [O IV]25.89/im luminosities la above 
those expected from their Lir. Likewise star- formation 
can explain the observed [O IV]25.89 /im luminosities for 
9 ±4% of the Seyfert galaxies (Figure [7]). Thus the 
[O IV] 25.89 /xm emission of Seyferts is generally domi- 
nated by the AGN, whereas for optically classified Hn 
galaxies star formation is likely to be the origin of the 
[O IV]25.89/im emission. 

5.2. The [Ne III] Excess 

The [Ne III] 15.56 /im emission may be produced by 
star formation and/or AGN. In this section we try 
to quantify the contribution of each one to the total 
[Ne III] 15.56 /im output. For this purpose we used the 
correlations between [Ne V]14.32 fim, [O IV]25.89/im and 
[Ne III] 15.56 /im (see Section 0]). Combining the typical 
[O IV]25.89/im/[Ne V]14.32/im and [Ne III] 15.56 /im/ 
[Ne V] 14.32 /im ratios found in AGN we estimate a ratio 
[Ne III]15.56/xm/[0 IV]25.89/im = 0.6, with rms scatter 
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Fig. 7. — [O IV]25.89 /im/Lia ratio vs. Lm luminosity. Galaxy 
symbols are as in Figure [5] The shaded region is the upper limit 
to the [O IV]25.89 /im/LjR ratio of star-forming galaxies. Galaxies 
above this region need an AGN contribution to the [O IV]25.89 fim 
emission in order to explain the observed ratio. 

of 0.3 for AGNs0, or log [Ne III]/[0 IV] = -0.22. This es- 
timated ratio is very similar to that obtained using a lin- 
ear fit between the [Ne III] 15.56 /im and [O IV]25.89/im 
luminosities (Table [3J . The median of this ratio for H II 
galaxies is 4.7, or log [Ne III]/[0 IV] = 0.67, which is 
~8 times larger than that estimated for active galax- 
ies. Figure [8] shows that the observed ratio for Seyfert 
galaxies and QSOs falls around the predicted value for 
AGN. Only a small fraction (4%) of the Hn galaxies 
have [Ne III] 15.56 fim/[0 IV]25.89/im ratios in the AGN 
range. It is also apparent that H n galaxies, a large frac- 
tion of LINERs and some Seyfert galaxies have an excess 
of [Ne III] 15.56 /im relative to their [O IV]25.89/im emis- 
sion which can also be attributed to star formation. 

We find that the four Hn galaxies with [Ne V] 14.32 /mi 
detections have the [Ne III]15.56 /im/[Ne V]14.32/im and 
[Ne II]12.81/im/[Ne V]14.32/im ratios much larger than 
those of Seyfert galaxies (Table [3|) . This indicates that 
the star-formation contribution is larger than that in 
Seyfert galaxies. Thus as suggested in Section[5]the AGN 
does not dominate the nuclear spectra of these galaxies. 

As can be seen in the left panel of Figure |U most of the 
outliers (galaxies more than 3cr above the [Ne V] 14.32 /im 
versus [Ne III] 15.56 /im flux correlation) have L[n c ml < 
10 42 erg s _1 . That is, for the most luminous ob- 
jects the AGN dominat es the [Ne II I ] 15. 56 /im emission. 
Using the relation by IHo fc Ketol (|2007f) that relates 
star formation rate (SFR) and the luminosity of the 
[Ne II] 12.81 /im and [Ne III] 15.56 /im emission lines and 
assuming a [Ne III]15.56/im/[Ne II]12.81/iin ratio = 0.3 
(see Section I5.3[) we find that the SFR of the outliers is 
between 0.1 and 100 M© yr _1 . 

5.3. The [Ne III] 15.56 iim/[Ne II]12.81 fim Ratio 

The [Ne III] 15.56 /im/[Ne II] 12.81 /im ratio traces the 
hardness of the radiation field an d the age of the stella r 
population (jRigbv fc Riekd 120041 : iSnijders et aTll2007h . 
Figure M compares the [Ne III]15.56/im/[Ne II]12.81 /im 
ratio observed in Hn galaxies (median 0.2) with that 

6 In this section the term "AGN" refers to galaxies classified as 
QSO or Seyfert. 
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JO IV]25.89 fim luminosity. Galaxy symbols are as in Figure 
[2 The solid red line is the [O IV]25.89 /im/[Ne IH]15.56/im 
ratio calculated from the linear fit to the [O IV]25.89 /im vs. 
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Fig. 9.— Distribution of the [Ne IIl]15.56^m/[Ne II]12. 81 /im ra- 
tio. The orange histogram corresponds to the observed ratios in H II 
galaxies, black corresponds to the observed ratios in active galaxies 
with star formation ([Ne III] 15.56 fim/[0 IV]25.89/mi ratios > 0.9). 
The blue histogram shows the ratios calculated for Seyfert galax- 
ies after subtracting the AGN contribution to the [Ne III] 15.56 (im 
emission. 

observed in Seyfert galaxies with star formation^ (me- 
dian 0.6). The larger ratio found in Seyfert galax- 
ies can be explained if the AGN contributes noticeably 
to the total [Ne III] 15.56 /an emission but not to the 
[Ne II] 12.81 /mi emission (see Section T4.2I) . By assum- 
ing this we calculate the star formation [Ne III] 15.56 /im/ 
[Ne II] 12.81 /im ratio in these Seyfert galaxies using the 
non-AGN [Ne III] 15.56 /mi emission that we estimated 
using the method described in Section [51 The median of 
this ratio becomes ^0.3, which is somewhat larger than 
that found in pure H II galaxies but it is i n the range of 
the o bserved ratio in this class of galaxies ([Brandl et al.l 
[20061) . 

6. COMPARISON OF STAR FORMATION IN SEYFERT 1 
AND 2 GALAXIES 

The AGN unification scenario predicts that there 
should be no differences between the star forma- 
tion activity in Seyfert 1 and 2 galaxies. A num- 
ber of studies compared the star formation rates in 

7 those with excess [Ne III] 15.56 ^m. 



Seyfert 1 and 2 galaxies. Some of them found en- 
hanced star formation activity in Seyfert 2 with re- 
spect Seyfert 1 galaxies, for a given AGN luminosity 
(IMaiolino et all 119951; IBuchanan et all 120061 IDeo et all 
l2007tlMelendez et al.rt2008bD , whereas others did not find 
any differences between the tw o types (jKauffmann et all 
120031; llmanishi fc Wadal I2004D . However, as discussed 
bv lShi et al.l ([2009D . the different methods used to study 
the star formation activity are sensitive to different stel- 
lar age ranges. Therefore these apparently contradic- 
tory results may be consistent with each other. The 
[Ne III] 15.56 /jm line is sensitive to young (<5Myr) mas- 
sive s tar formation ()Rigbv fe Riekd [20041 ; ISniiders et all 
120071 ). Studies using this line (and the [Ne II] 12.81 /mi 
line) found t hat star formati o n is enhanced in Seyfer t 
2 galaxies (IDeo et all [20071; IMelendez efaLl I2008bl) . 
Kauff mann et all (|2003l) did not find differences be- 
tween Seyfert 1 and Seyfert 2 star formation, but their 
method (the 4000A break and the H<5 stellar absorption) 
was sensitive to stellar populations older than 100 Myr 
(IShi et alll2009D. The same appl ies to the PAH features 
used bv llmanishi fc Wadal {2004), which can be affected 
by th e presence of the AGN (|Diamond-Stanic fc Riekel 
2010, and references therein). 

It is also important to use the appropriate indicator 
for the AGN luminosity. We will use the [O IV]25.89 (im 
emission, as it is a good isotropic indicator (jRigbv et all 
I2009t iDiamond-Stanic et al J 120091) . as long as the AGN 
dominates the [O IV] 25. 89 /im emission as we showed in 
Section 15. II 

We estimate the star formation contribution to 
the [Ne III] 15.56 /mi emission by using the observed 
[OlV]25.89/zm/[NeIII]15.56/mi ratio and subtract- 
ing the estimated AGN contribution from the total 
[Ne III] 15.56 /im line strength. Further details 



are 



O IV 
O IV 



25.89 /mi/ 
25.89 /mi/ 



provided in Appendix [AJ We used the 
[Ne III] 15.56 /im instead of the 
[Ne II] 12.81 /im ratio because we could not calcu- 
late the typical AGN [O IV]25.89/mi/[Ne II]12.81 /im 
ratio due to the larger relative contribution of star 
formation activity to the [Ne II] 12.81 /mi emission 
(Section I4.2[) . We only subtract this contribution from 
those active galaxies more than la above the typical 
AGN [Ne III]15.56/im/[0 IV]25.89/mi ratio (see Section 
15.21 and Figure [5]) , since the galaxies below this ratio are 
presumably dominated by the AGN and their star forma- 
tion activity, if any, is masked. This criterion basically 
selects only those active galaxies with at least 25% of the 
[Ne III] 15.56 /im emission arising from star formation. 
In our sample 45% of the AGNs with the [O IV]25.89 /mi 
line detected have excess [Ne III] 15.56 /im emission due 
to star formation. The right panel of Figure [TU] shows 
that the fraction of Seyfert 2s with star formation is 
slightly higher than that of Seyfert Is. In all the lumi- 
nosity bins except for one, the difference in the fraction 
with a [Ne III] 15.56 /im excess is less than 2a significant. 
However, for the overall proportion of [Ne III] 15.56 /im 
excess detections, the Fisher's exact test indicates that 
there is a probability of less than 0.03 that the incidence 
of this behavior is the same in the Syl and Sy2 samples. 
Therefore, the full sample in our study supports, with 
a moderate statistical significance, previous indications 
that Sy2 host galaxies tend to have higher rates of star 



9 



SFR detection limit (M yr 1 ) 
0.1 1 10 100 1000 



20 
10 




Sy1 




38 39 40 41 42 43 
log L[OIV] (erg s" 1 ) 




[] 



Sy1 a 
Sy2 □ 



4± 



38 39 40 41 42 43 
log L[OIV] (erg s" 1 ) 



38 39 40 41 42 43 
log L[OIV] (erg s" 1 ) 



Fig. 10. — Distribution of the [O IV]25.89/im luminosities of the Seyfert 1 (left) and 2 (middle) galaxies. The filled histograms are the 
galaxies with a [Ne III] 15.56 fim excess. The right panel compares the fraction of type 1 and type 2 Seyferts with a [Ne III] 15.56 /im excess 
for each luminosity bin. We used Poisson statistical errors to estimate the uncertainties of the fractions. The scale on the upper part of 
the right panel gives the lower limit to the SFR that can be detected with our method for each AGN luminosity bin. 



formation than Syl hosts. However the full sample is 
heterogeneous and may be subject to a number of biases. 
For example, it includes a number of radio galaxies, 
which are known to have low star formation rates - in 
fact none of the radio galaxies have [Ne III] 15.56 /mi 
excess. It also contains ULIRGs with Seyfert spectral 
characteristics, most of which do have [Ne III] 15.56 /im 
excess. The observed differences depend to some extent 
on the relative number of galaxies in these two classes. 
In the following section, we discuss this trend in two 
relatively complete but smaller subsamples where the 
biases should be reduced. 

6.1. Results for the 12jim sample and the RSA samples 

At this point we study separately the 12 /im sample and 
the RSA sample. Both are complete samples of Seyfert 
galaxies selected with homogeneous criteria. The former 
is selected based on their IRAS 12 /jm fluxes whereas 
the latter is selected based on the optical magnitude of 
the host galaxy. Given the different selection criteria 
they they might be affected by star formation in different 
ways. 

As can be seen from Table 0] and Figure [TTJ the frac- 
tion of galaxies with a [Ne III] 15.56 /tm excess in the 
RSA sample is higher (by a factor of 1.4) than in the 
12 /im sample. The lower number of galaxies with a 
[Ne III] 15.56 /tm excess in the 12 /im sample is because 
the galaxies of this sample are brighter and the sensitiv- 
ity to star formation of our method is reduced at high 
[O IV]25.89/im luminosities. Both effects are clearly seen 
in Figure [TTJ First, the RSA sample contains a larger 
fraction of low luminosity Seyferts (L[ IV ] < 10 40 erg 
s _1 ) than the 12 /im sample. Second, for each AGN lu- 
minosity, we are sensitive to a different SFR limit. In 
other words, for the brightest Seyferts we are only sensi- 
tive to SFR above 100 yr" 1 , while for low luminosity 
AGN we are sensitive to SFR > 0.1 M Q yr" 1 . We note 
that the fraction of galaxies with star formation in each 
luminosity bin is similar for both samples (right panel of 
Figure ITT)) . 

In both samples the fraction of Seyfert 2s with 
[Ne III] 15.56 /tm excess is slightly larger than that of 
Seyfert Is. However due to the small size of the samples 
the statistical significance of this excess is low. Based on 
the two sample K-S test, the probability for these differ- 
ences being due to chance is 0.25 and 0.4 for the 12 /tm 



TABLE 4 
Fraction of galaxies with a 
[ne iii] excess in the rsa and 

12 flM SAMPLES 



Type 


12 /tm sample 


RSA 


Seyfert 1 


15/45 


12/25 


Seyfert 2 


23/52 


25/45 


Total 


38/97 


37/70 



Note. — In this ta- 

ble we only include galaxies 
with the [O IV]25.89/mi and 
[Ne III] 15.56 /tm emission lines de- 
tected. Some galaxies are common 
to both samples, 
and RSA samples respectively. 

7. LINE RATIOS AND MODELS 

We studied the physical conditions in the narrow line 
region of t he AGNs with the p hotoionization code MAP- 
PINGSIII ([Groves et alJ2004) . We used a radiation pres- 
sure dominated model that includes the effects of dust. 
For the i nput paramet e rs we followed the prescription 
given by iGroves et all (|2006l ). Briefly, we assumed a 
plane parallel geometry and solar abundances. We mod- 
eled the input ionizing sp ectrum with two p ower-laws 
with exponential cut-offs (jNagao et al.l [2001). We ex- 
plored the effect of the variations in the total hydrogen 
column density ranging from log nn (cm -2 ) = 19.0 to 
22.0. This range was chosen since it reproduces the ob- 
served line ratios and is in good agreement with the hy- 
drogen colum n density determined using UV and X-ray 
observations ([Crenshaw et al.l 120031) . For each hydrogen 
column density we varied the total pressure (P/k) and 
the incident ionizing flux (Jo)- The values for the total 
pressure are log P/k (K cm~ 3 ) = 6, 7, 8, 9. These cor- 
respond to electron densities, as traced by [S IL]A6716A/ 
A6731A, of - <10\ 10 2 , 10 3 and > 10 4 cm 3 respec- 
tively. The explored range for the ionizing flux is Iq — 
0.1 to 0.55. The incident ionizing flux is scaled by the fac- 
tor 2.416 x((P/fc)/10 6 ) erg cm~ 2 s _1 which gives a range 
in the ionization parameter, log U, from ~ -3 to -2. 

7.1. The [Ne V] ratio 
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Fig. 11. — Di strib ution of the [O IV]25.S9 fim luminosities for the RSA sample (left panel) and 12 fim sample (middle panel). Symbols 
are as in Figure ITul 



Figure [T2] shows the ratio between the two [Ne V] 
emission lines available in the mid-infrared. The criti- 
cal densities of the [Ne V] 14.32 /jm and [Ne V]24.32^im 
lines are 3.5xl0 4 cm~ 3 and 6.2xl0 3 cm~ 3 , respectively 
(jOsterbrock fe Ferlandll200"6lh thus they can be used to 
trace the density of the NLR. We find however that a 
large fraction of the galaxies (~ 30%) lies above the low 
density limit ratio. Th is issue and i t s pos sible cause are 
discussed in detail by IDudik et al.l (|2007l ) and also ad- 
dressed b ylTommasin et al.l (12010ft for a sample of Seyfert 
galaxies. Dudik et al.l ([20071 ) concluded that the ratios 
above the low density limit are due to differential extinc- 
tion from the obscuring torus, and thus can indicate the 
in clination ang l e to o ur line of sight. However, similar 
to lBaum et aTJ(|2010l ). we do not find any significant dif- 
ference between the proportion of Seyfert 1 and Seyfert 
2 galaxies above the low density limit, as would be pre- 
dicted by this hypothesis. 

We calculated the A24.3/A14.3 ratio to quantify the ex- 
tinction needed to move the points above the low den- 
sity limit ratio to this value (which would be a lower 
limit to the extinction since the real ratio may be lower 
than the low density limit). The ratios are summarized 
in Table [5] for some infrared extinction laws. As can 
be seen depending on the extinction law chosen the ex- 
tinction can increase or decrease the [Ne V]24.32/im/ 
[Ne V] 14.32 fim ratio. Since the construction of an ex- 
tinction law implies some interpolation we also used di- 
rect m easurements of A 2 4 / Ak and Ai 5 / Ak Uiang et ah! 
(2006) measured the extinction at 15 //m and found A\^j 
A K = 0.40, although it ranges from 0.25 to 0.55. For 
the A24/ A.K ratio the range is 0.28 to 0.65 for different 
Ak bins (Cha pman et al.ll2009D and the average value is 
~0.5 ([Chapman et al.ll2009t iFlahertv et al.ll2007f ). From 
these measurements we estimate A24.3/A14.3 — 1.2 ± 0.3. 
Again it is close to unity and it is not clear if the extinc- 
tion would increase or decrease the [Ne V] ratio. More- 
over the extinction is nearly neutral between the lines 
with all the extinction laws and therefore there is very 
little potential effect on the ratio unless the extinction is 
extremely large. 

We also checked if aperture effects can explain the 
large values of this ratio. The [Ne V] 14.32 /im line 
is observed with the SH module (5'.'7x 11'.'3), whereas 

8 From the extinction laws, the difference between these A values 
and those at 14.3 fim and 24.3 fim would be small (~5%). 



TABLE 5 

Extinction at 24 fim relative 

TO THAT AT 14 fim 



Extinction law 


A24/A14 


Rosenthal" 1 


1.25 


McClure 0.3 < A K < l b 


0.83 


McClure 1 < A K < 7 b 


0.83 


Chiar Galactic Center 


1.38 


Chiar Local ISM C 


0.97 



b McClurdl200a 



Chiar fe Tieiensl [200(3 

the [NeV]24.32//m is observed with the LH module 
(ll / .'lx22'.'3). We would expect a correlation between the 
[Ne V]24.32^m/[Ne V]14.32^m ratio and the distance if 
the [Ne V] emission region was larger than the SH slit. 
However we do not find this correlation thus we rule out 
any aperture bias. 

Alternatively the atomic parameters for the Ne + ion 
may not be sufficiently accurate for the comparison with 
these observations and the calculated low density limit 
may be incorrect. 

7.2. Ionization parameter 

In Figure [12] we plot the model grid together with the 
data. In this plot we vary the ionizing flux, irji and the 
total pressure, P/k, for a constant total hydrogen column 
density (log n H (cm" 2 ) = 20.3). We see that the [Ne V] 
ratio can be used to measure the pressure. However the 
pressure (or density) values derived from this ratio have 
some issues, since a considerable fraction of galaxies are 
above the low density limit ratio (Section [73}. 

Also, Figure [H shows that the [O IV]25.89 / um/ 
[NeV]24.32^m ratio traces the ionization parameter. 
The [O IV]25.89/im/[Ne V]24.32^m ratio spans a nar- 
row range (less than 1 dex) which is related to the 
tight correlation found between the [O IV]25.89/im and 
[Ne V]24.32 /im luminosities (Section [4]). This implies 
a very small range for the ionization parameter (-2.8 
< logU < -2.5). This range is similar to that found 
when modeling the o ptical [O HI] and H /? emission 
(|Baskin fc Laorj 120051: Kraemer et atl I1999D . The crit- 
ical density of the [O III] lines is larg er than that of 
these mid-IR lines and for instance the iBaskin fc Laorl 
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Fig. 12.— Predicted [Ne V]24.32 /im/[Ne V] 14.32 ^tm ratio vs. 
[O IV]25.89/im/[Ne V]24.32/im. Galaxy symbols are as in Figure 
[2] The dashed black line indicates the low density limit calculated 
at 10000 K for the [Ne V]24.32 /mi/[Ne V]14.32 /im ratio (galaxies 
above this line have densities below the low density limit). The 
blue lines trace constant pressure and the orange lines constant 
ionization parameter (Jo). The input parameters for the model are 
log P/k = 6, 7, 8, 9 and I = 0.1, 0.15, 0.20, 0.25, 0.32, 0.4, 0.55 
for constant hydrogen column density log ran (cm -2 ) = 20.3. The 
galaxy symbols are as in Figure [2] 

(2005) models predict electron densities logn e (cm -3 ) = 
5.8 ±0.7 well above the critical density of the mid-IR 
lines studied here. Thus it is possible that these lines are 
produced in different gas clouds. 

7.3. Hydrogen column density 

We use the [Ne III] 15.56 /xm and the [Ne V]14.32/mi 
lines to estimate the total hydrogen (H I + H II) 
column density, njr. The [Ne V] 14.32 /im line is 
produced in the inner part of the NLR while the 
[Ne III] 15.56 /im is produced in a more external region. 
Thus the [Ne III]15.56/mi/[Ne V]14.32/mi ratio is sub- 
ject to change depending on the column density. As 
expected, the models predict larger [Ne III] 15.56 /im/ 
[Ne V] 14.32 /im ratios for larger n^. 

In Figure [13] we plot the model grid together with the 
observed ratios. The [Ne III]15.56/im/[Ne V]14.32/xm 
ratio is not sensitive to column densities larger than log 
71h (cm -2 ) > 21. We find that the column density for 
most of the galaxies is in the range 20.3 < log Tin (cm -2 ) 
< 21. There is no difference between the predicted hy- 
drogen column density for type 1 and 2 galaxies. This 
column density range is compatible with the values de- 
rived from UV observation s although lower than that de- 
termined using X-ray data (|Crenshaw et al.ll2003|) . A few 
galaxies are above the model grid. The [Ne III] 15.56 /im/ 
[Ne V] 14.32 /im ratio of these galaxies is larger than that 
predicted by the models. As we proposed in Sections 14.21 
and 15. 2\ star formation acting in these galaxies is likely 
producing the extra [Ne III] 15.56 /im emission. 

8. CONCLUSIONS 

We studied a sample of 426 galaxies observed 
with the Spitzer /IRS spectrograph in the high res- 
olution mode. Our analysis includes published 
data for QSO, Seyfert, LINER and Hn galaxies 
as well as unpublished measurements for the RSA 
Seyfert sample and LIRGs. We explored the rela- 
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Fig. 13.— Predicted [Ne IH]15.56/mi/[Ne V] 14.32 fim ratio vs. 
[O IV]25.89/im/[Ne V]24.32 fim. The green lines trace constant hy- 
drogen column density and the orange lines constant ionization pa- 
rameter. The input parameters for the model are log ran (cm - 2 ) 
= 19.0, 20.0, 20.3, 20.6, 21.0, 22.0 and I = 0.1, 0.15, 0.20, 0.25, 
0.32, 0.4, 0.55 for constant pressure (log P/k = 8). The galaxy 
symbols are as in Figure [2] 

tionship between the high-ionization ([Ne V] 14.32 /im, 
[Ne V]24.32/im and [O IV] 25.89 jum) and intermediate 
ionization ([Ne III] 15.56 /im and [Ne II] 12.81 /mi) emis- 
sion lines present in the mid-IR spectra. The median 
ratios for each class of galaxies are listed in Table [3] The 
main results are as follows: 

1. There is a tight linear correlation between the 
fluxes and luminosities of the high-ionization emis- 
sion lines [O IV]25.89 fim and [Ne V]24.32 /mi. This 
correlation spans 5 orders of magnitude (38 < 
log i[Ncv]24(erg s -1 ) < 43), ranging from LIN- 
ERs, low-luminosity Seyfert galaxies, types 1 
and 2, to QSOs. The typical [O IV]25.89/mi/ 
[Ne V]24.32/im ratio for AGNs is 3.5, with rms 
scatter of 0.8. The correlation also holds be- 
tween the [O IV]25.89/im and the [Ne V]14.32/mi 
lines, although the scatter is larger. The typical 
[O IV]25.89/mi/[Ne V]14.32/im ratio for AGNs is 
3.4, with rms scatter of 1.4. 

2. There is also a good linear correlation between 
the [Ne III] 15.56 /jm and the [Ne V] 14.32 //m fluxes 
and luminosities for Seyfert galaxies and QSO. 
The [NeIII]15.56/zm/[Ne V] 14.32 /zm ratio is 1.9, 
with rms scatter of 0.8. Using this correlation 
and the previous one we calculated the typical 
[Ne III]15.56/zm/[0 IV]25.89/im ratio for AGNs of 
0.6, with rms scatter of 0.3. 

3. We calculate the [O IV]25.89/im emission due 
to star formation (logLp i V ] (erg s -1 ) < 
logLiR (erg s _1 ) — 4.7 ± 0.6) and we estimate 
that it may dominate the total [O IV]25.89/mi 
emission when the intrinsic AGN luminosity is a 
factor of 20 smaller than the star forming lumi- 
nosity. In general, we find that the [O IV]25.89/im 
emission is dominated by the A&N in Seyfert 
galaxies, whereas star-formation can explain the 
[O IV]25.89/im emission of optically classified Hn 
galaxies. 
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We do not find any significant difference between 
the mid-IR high-ionization emission lines in type 
1 and in type 2 Seyfert galaxies. Either there are 
no differences in the conditions in the NLR of the 
two AGN types, or the effects associated with the 
different line of sight (i.e., dust extinction) are min- 
imized in the mid-infrared spectra. 

We find that a significant number of Seyfert galax- 
ies from the 12 /im sample (30%) and the RSA sam- 
ple (40%) show evidence for excess [Ne III] 15.56 /im 
emission relative to their [O IV] 25. 89 /im emission 
associated with star formation. The larger fraction 
in the RSA sample is explained because the RSA 
sample contains a larger fraction of low-luminosity 
AGNs (L[q iv i < 10 41 erg s _1 ) in which these star 
formation excesses are easier to detect. 



The fraction of Seyfert 2 galaxies with 
[Ne III] 15.56 /xm excess is larger than that of 
Seyfert 1 with a moderate statistical signicance 
in the full sample although this sample may be 
affected by several bias. However when we consider 
the 12 /im or RSA Seyfert samples separately these 
differences are not statistically significant due to 
the smaller size of these samples. 

A considerable fraction (30%) of the galaxies have 
[Ne V]24.32/im/[Ne V] 14.32 /xm ratios above the 
low density limit. We did not find a connection 
between the Seyfert type and the ratio of galaxies 
above the limit. We are not able to explain this in 
terms of differential extinction. 

Our modeling shows that the nebular conditions 
in the NLRs are remarkably similar among all the 



AGNs in our sample. This similarity allows us to 
compare conditions critically in the NLRs of the 
type 1 and type 2 galaxies. There appear to be no 
significant overall differences, consistent with the 
unified model. We constrained the ionization pa- 
rameter in the range -2.8 < log U < -2.5 and the 
hydrogen column density 20 < log fin (cm -2 ) < 21. 

The relationships presented in this paper provide an 
important benchmark for the interpretation of the future 
mid-IR observations of AGNs and star forming galaxies 
with JWST/MIRI. 
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APPENDIX 
CALCULATING THE AGN CONTRIBUTION 



In this appendix we explain briefly the method used in Sections [5] and [6] to calculate the star formation fraction of 
the [Ne III]15.56/im emission. 

The total [Ne III] 15.56 /xm emission includes two components, one from the AGN and other from star formation 
(SF). 

[Ne in] = [Ne iii] sf + [Ne iii]agn (Al) 

If we assume that all the [O IV]25.89/im is produced by the AGN we can use the typical [Ne III] 15.56 /im/ 
OlV]25.89/im ratio observed in Seyfert galaxies with low star formation (Figure [8]) to estimate the amount of 
Ne III]15.56/im emission coming from the AGN. 



[Ne in] = [Ne iii] S f + [O iv] 



[Ne in] 
lOrTf 



AGN 



Finally we obtain the fraction of [Ne III] 15.56 /xm from star formation. 

[Ne in] SF 1 [O iv] /[Nem] 



[Ne in] 



1 



[Ne in] V [O iv] 



(A2) 



(A3) 



AGN 



This method can be applied as well to other line ratios. The only assumption is that one of the lines in the 
ratio is uniquely produced by the AGN while star formation and the AGN contributes to the other line (e.g. the 
[Ne III]15.56/im and [Ne V] lines). Also note that this method makes use of a typical AGN ratio thus the estimated 
star formation contribution to an emission line for a single object will be uncertain and dependent on the ionization 
paremeter. 
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9.27 


Svl 

oy -L 


17 


11 


16 


10 


32 


2 


TRAS F120184-1941 

X X V ii O J_ X ,i X O 1^ X J^t X 


12 


04 


24 


.5 


+ 19 


25 


09 


813 


LINER 


3.0 


<0.80 


0.35 


<0.10 


<0.63 


3 


NGC 4085 


12 


05 


22 


.7 


+50 


21 


10 


14.6 


H 11 


23 


<0.15 


2.9 


<0.92 


<0.44 


7 


NGC 4088 


12 


05 


34 


.2 


+50 


32 


20 


13.4 


H 11 


37 


<0.39 


2.5 


<0.50 


0.73 


8 


NGC 4125 


12 


08 


06 


.0 


+65 


10 


26 


22.9 


LINER 


2.3 


<0.62 


3.3 




0.75 


5 


IRAS 12071-0444 


12 


09 


45 


.1 


-05 


01 


13 


600 


Sy2 


5.2 


2.9 


5.1 


3.7 


6.6 


3 


NGC 4151 c d 


12 


10 


32 


.6 


+39 


24 


20 


14.0 


Svl 

oy 1 


131 


77 


205 


68 


244 


2 


NOn 41 57 


12 


11 


04 


.4 


+50 


29 


04 


13.3 


H 11 


12 


<0.070 


1.4 


<0.70 


1.1 


7 


NGC 4192 


12 


13 


48 


.3 


+ 14 


54 


01 


19.3 


H 11 


19 


<0.18 


4.7 


<0.48 


1.6 


2 


NGC 4194 


12 


14 


09 


.5 


+54 


31 


36 


40.3 


Sy2 


165 


3.0 


54 


4.0 


27 


8 


PG 1211 + 143 


12 


14 


17 


.7 


+ 14 


03 


12 


368 


QSO 


0.32 


<0.43 


0.73 


<0.56 


2.4 


1 


NGC 4236 


12 


16 


42 


.1 


+69 


27 


45 


4.36 






<1.5 






<0.23 


5 


NGC 4235 c 


12 


17 


09 


.9 


+07 


11 


29 


33.5 


Syl 


3.8 


<0.43 


3.3 


<0.54 


3.3 


2 


Mrk 766 d 


12 


18 


26 


5 


4-29 


48 


46 


55.9 


Svl 
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23 


21 


24 


18 


46 


4 
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12 


18 
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24 


59 
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5 
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12 


18 
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.5 


+47 


18 


14 
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12 


<1.3 
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2 
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12 


20 
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16 
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3.9 


<0.29 


0.88 


2 
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12 


22 


54 


.9 
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49 


20 
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5 
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12 


25 


42 


.8 


+00 


34 


21 


30.1 


H 11 


29 


<0.57 


3.8 


<0.83 


1.3 


8 


NGC 4388 c d 


12 


25 


46 


.7 


+ 12 


39 


43 


36.3 


Sy2 


75 


45 


106 


68 


308 


2 


NGC 4395 c 


12 


25 


48 


.9 


+33 


32 


48 


4.33 


Syl 


4.9 


0.93 


6.8 


1.4 


6.9 


2 


NGC 4419 


12 


26 


56 


.4 


+ 15 


02 


50 


17.2 


H 11 


36 


<0.46 


8.5 


<1.0 


<1.6 


2 
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12 


27 


40 


.5 
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04 


44 
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H 11 


6.8 


<0.24 


2.1 
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2 


NGC 4450 


12 


28 


29 


.6 
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05 


05 


16.6 


LINER 


3.5 


<0.53 
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0.63 


5 


NGC 4457 


12 


28 


59 


.0 


+03 


34 


14 


12.4 


LINER 


7.8 


<0.18 


3.2 


<0.47 


2.2 


2 
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12 


29 


06 


.7 


+02 


03 


08 
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Syl 
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3.4 


6.0 


2.9 
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3 


NGC 4477 c 


12 


30 
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38 


11 
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3.2 
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TABLE 1 — Continued 



Name 


R.A. a 
(J2000.0) 


Decl. a 
(J2000.0) 


Dist. 

(Mpc) 


Spcct. 
Class. a 


[Ne ii] 
12.81 fim 


[Ne v] 
14.32 /im 


[Ne in] 
15.56 /im 


[Nc v] 
24.32 /im 


[O iv] 
25.89 /im 


Ref. 
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7 
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12 


30 


49 


.4 
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23 


28 
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1.2 


2 
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12 


31 


59 


2 
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25 


13 
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2 
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12 


32 


03 


6 
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09 


29 
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0.91 


1.3 
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1 


NGC 4536 


12 


34 


27 


.1 
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11 


17 
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<7.9 


5 


NGC 4507 c 


12 


35 


36 


.6 
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54 


33 


50.9 


Sy2 


30 


12 


29 


8.6 


34 


2 


NGC 4552 


12 


35 


39 


.8 
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33 


22 


15.9 


LINER 


1.7 


<0.89 
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<0.53 


5 


NGC 4559 


12 


35 


57 


.7 


+27 


57 


35 
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H ii 
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0.12 


5 


NGC 4565 c d 


12 


36 


20 


.8 
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59 


15 
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Syl 


2.6 


0.32 


3.3 
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2 


NGC 4569 
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36 


49 
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09 


46 
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36 
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5 
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12 


37 
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05 
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Syl 
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<0.40 


12 
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2 
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12 
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20 
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2 


NGC 4594 c d 


12 
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Syl 


14 
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12 
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07 


58 
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Syl 
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0.82 
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<2.3 


4 
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12 


40 


52 


.9 
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45 


21 
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42 
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3.5 


17 


2 


NGC 4625 


12 


41 


52 


.7 
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16 


26 


9.20 
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<0.84 
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<0.26 


5 


NGC 4631 


12 


42 
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29 


8.10 
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<0.52 
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5 
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12 


42 


52 


.4 
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15 


26 
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Syl 


0.84 


<0.27 


1.1 


<0.28 
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2 


NGC 4666 


12 


45 


08 
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27 


42 
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LINER 


35 
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8.3 


1.4 
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2 


NGC 4676 


12 


46 


10 


.1 


+30 


43 


55 


94.0 


H ii 


28 


<0.25 


4.7 


<0.33 


1.4 


8 


PG 1244+026 


12 


46 


35 


2 


+02 


22 


08 


213 


QSO 


0.94 


0.53 


1.2 


1.1 


1.5 


1 


NGC 4725 c d 


12 


50 


26 


6 


+25 


30 


03 


17.4 


Sy2 


1.2 


<0.38 


2.4 


<0.50 


1.8 


2 


NGC 4736 


12 


50 


53 


.1 


+41 


07 


13 


5.00 


LINER 


13 


<0.81 


14 




3.3 


5 


NGC 4748 d 


12 


52 


12 


.5 


-13 


24 


53 


63.4 


Syl 


7.4 


6.7 


16 


20 


82 


4 


IRAS F12514+1027 


12 


54 


00 


.8 


+ 10 


11 


12 


1667 


Sy2 


2.3 


1.9 


2.7 


1.7 


2.7 


3 


Mrk 231 


12 


56 


14 


.2 


+ 56 


52 


25 


186 


Syl 


20 


<3.0 


3.0 


<18 


<9.5 


3 


NGC 4826 


12 


56 


43 


.7 


+21 


40 


57 


5.00 


Sy2 


105 


<1.3 


23 




4.3 


5 


NGC 4818 


12 


56 


48 


.9 


-08 


31 


31 


9.40 


H ii 


185 


<1.2 


14 


<2.0 


<1.6 


8 


NGC 4922 d 


13 


01 


24 


.9 


+29 


18 


40 


103 


LINER 


36 


2.4 


9.2 


<1.9 


4.3 


2 


MCG -02-33-098-E 


13 


02 


19 


.7 


-15 


40 


03 


68.9 


H ii 


32 


<0.58 


12 


<0.61 


1.2 


2 


MCG -02-33-098-W 


13 


02 


20 


.4 


-15 


45 


59 


68.1 


H ii 


67 


<0.51 


12 


<0.79 


<1.2 


2 


NGC 4941 <= d 


13 


04 


13 


1 


-05 


33 


05 


15.9 


Sy2 


13 


7.2 


24 


6.9 


26 


2 


NGC 4939 c 


13 


04 


14 


.4 


-10 


20 


22 


44.8 


Sy2 


8.1 


12 


24 


17 


65 


2 


NGC 4945 c 


13 


05 


27 


.5 


-49 


28 


05 


8.19 


Sy2 


586 


3.4 


69 


<5.1 


35 


2 


PG 1302-102 


13 


05 


33 





-10 


33 


19 


1422 


QSO 


0.36 


0.49 


0.66 


0.39 


2.6 


1 


UGC 8201 


13 


06 


24 


.9 


+67 


42 


25 


3.68 






<0.75 






<0.24 


5 


NGC 4968 d 


13 


07 


06 





-23 


40 


37 


42.6 


Sy2 


25 


18 


34 


11 


34 


4 


PG 1307+085 


13 


09 


47 


.0 


+08 


19 


48 


739 


QSO 


0.40 


0.56 


0.98 


0.62 


0.74 


1 


NGC 5005 c d 


13 


10 


56 


2 


+37 


03 


33 


13.8 


Sy2 


37 


<1.1 


13 


<1.2 


4.6 


2 


PG 1309+355 


13 


12 


17 


.8 


+35 


15 


21 


893 


QSO 


0.51 


0.27 


1.3 


<0.30 


<0.49 


1 


NGC 5033 c d 


13 


13 


27 


.5 


+36 


35 


38 


12.4 


Syl 


31 


1.2 


14 


2.0 


9.2 


2 


IC 860 


13 


15 


03 


.5 


+24 


37 


07 


55.8 


H ii 


3.6 


<0.14 


0.68 


<0.48 


<0.62 


2 


IRAS 13120-5453 


13 


15 


06 


.3 


-55 


09 


22 


136 


Sy2 


150 


1.7 


18 


<20 


6.4 


3 


NGC 5055 


13 


15 


49 


3 


+42 


01 


45 


7.80 


LINER 


22 


<0.95 


9.6 




2.3 


5 


UGC 8387 


13 


20 


35 


.3 


+34 


08 


22 


101 


H ii 


113 


<0.54 


19 


<2.3 


6.2 


2 


NGC 5104 


13 


21 


23 


.1 


+00 


20 


32 


80.0 


LINER 


43 


<0.39 


5.1 


<1.2 


2.4 


2 


MCG -03-34-064 d 


13 


22 


24 


.5 


-16 


43 


42 


71.7 


Syl 


56 


63 


119 


38 


115 


4 


NGC 5128 c 


13 


25 


27 


.6 


-43 


01 


08 


7.85 


Sy2 


197 


23 


150 


28 


135 


2 


NGC 5135 c d 


13 


25 


44 


.1 


-29 


50 


01 


58.9 


Sy2 


112 


14 


58 


18 


73 


2 


NGC 5194 c d 


13 


29 


52 


.7 


+47 


11 


42 


6.77 


Sy2 


44 


3.0 


34 


3.9 


19 


2 


NGC 5195 


13 


29 


59 


.6 


+47 


15 


58 


8.00 


LINER 


18 


<1.0 


6.9 




<1.9 


5 


3C 286 


13 


31 


08 


.3 


+30 


30 


32 


5399 


Syl 




0.75 


0.69 






6 


NGC 5218 


13 


32 


10 


.4 


+62 


46 


03 


41.1 


H ii 


45 


<0.16 


4.9 


<0.57 


<2.3 


2 


MCG -06-30-015 d 


13 


35 


53 


.8 


-34 


17 


44 


33.4 


Syl 


5.0 


5.0 


5.9 


7.4 


26 


4 


IRAS F13342+3932 


13 


36 


24 


.1 


+39 


17 


31 


866 


Syl 


5.7 


3.5 


5.0 


4.2 


10 


3 


NGC 5236 


13 


37 


00 


.9 


-29 


51 


55 


3.60 


H ii 


503 


<0.61 


29 


<1.2 


5.7 


7 


Mrk 266 d 


13 


38 


17 


.5 


+48 


16 


37 


116 


Sy2 


57 


8.0 


28 


11 


53 


8 


NGC 5273 c 


13 


42 


08 


.3 


+35 


39 


15 


15.2 


Syl 


2.0 


<1.5 


3.4 


<0.88 


4.9 


2 


Mrk 273 


13 


44 


42 


.1 


+55 


53 


12 


167 


LINER 


42 


12 


34 


15 


56 


3 


IRAS F13451 + 1232 


13 


47 


33 


.4 


+ 12 


17 


24 


565 


Sy2 


5.0 


<1.0 


5.1 


<2.1 


2.1 


3 


IC 4329A d 


13 


49 


19 


.3 


—30 


18 


34 


69.6 
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28 


29 


57 


35 


117 


4 


PG 1351+640 


13 


53 


15 


.8 


+63 


45 


45 


402 
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1.8 


<0.92 


2.7 


<1.00 
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1 
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13 


53 
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29 


26 
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13 


55 


39 
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27 


42 


37 1 
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2 
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13 


56 


02 


.9 
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22 


19 
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Svl 
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9.3 


18 


41 


20 


69 


3 


NGC 5395 c 


13 


58 


38 





+37 


25 


28 


50.3 


Sy2 




<0.54 




<0.31 


<0.18 


2 


NGC 5398 


14 


01 


21 


.6 


-33 


03 


49 


10.8 






<0.50 






<0.31 


5 


NGC 5408 


14 


03 


20 


.9 


-41 


22 


39 


4.85 


H ii 




<0.93 


1.2 




<0.51 


5 


NGC 5427 c 


14 


03 


26 


.1 


-06 


01 


50 


37.5 


Sy2 


10 


1.7 


5.1 


1.2 


4.2 


2 


NGC 5474 


14 


05 


01 


.6 


+53 


39 


44 


6.40 


H ii 




<0.83 






<0.31 


5 


IRAS F14070+0525 


14 


09 


31 


3 


+05 


11 


31 


1341 


Sy2 


1.8 


<0.15 


0.58 


<0.60 


<1.5 


3 


3C 295 


14 


11 


20 


.6 


+52 


12 


09 
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Sy2 


0.060 


0.070 
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6 


Circinus c 


14 


13 


09 


.9 


-65 


20 


20 


6.09 


Sy2 
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2 


NGC 5506 c d 


14 


13 
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12 
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81 


58 
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63 
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2 


PG 1411+442 


14 


13 
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00 


13 
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0.96 


0.92 


0.55 


1.5 


1 


NGC 5548 d 


14 


17 


59 


.5 
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08 


12 


74.5 


Syl 


8.5 


5.4 


7.3 


3.9 


17 
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R.A. a 
(.T2000.0) 
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14 


45 


09 


1 


-20 


52 


13 
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H 11 


77 
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8.2 


<0.23 


0.93 


2 


NGC 5743 


14 


45 
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-20 


54 


48 
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H 11 


50 
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0.78 


3.2 


2 
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14 


51 
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09 


26 
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0.51 
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1 
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14 


57 


41. 
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07 
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25 
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24 
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14 


57 


44. 


6 


-43 


07 


53 


65.4 
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18 
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14 
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6 
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40 


19 
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6 
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References . — (1 ) IVeilleux et al.1 120091). (2) This w ork. f31IFarrah et alJ 120071). < 4 ) [Tommasin et al.l 120081120101 ). f5) IDale et~aT1 120091) . (6) 
IGoriian et al.l 120071) . (7) [Goulding fc Alexanderf 120091) . (8) IBernard-Salas et al.1 120091) 

Note. — Fluxes are expressed in units of 10 -14 erg cm -2 s — . 

a Coordinates and optical spectroscopic classification from NED. 

^ Wc calculated the distance from the redshift assuming a flat cosmology with Hq — 70 km s _1 Mpc _1 , Om — 0.3, and Oa — 0.7. Except for the 
galaxies from Dale ct al. (2009), Bcrnard-Salas et al. (2009) and Goulding & Alexander (2009) for which wc used the distances adopted by these 
authors. 

c Member of the RSA sample. 
^ Member of the 12 jum sample. 
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